PROCEEDINGS 

OF SCIENCE 



Delta 1=3/2 K to pi-pi decays with nearly physical 
kinematics 



o 

(N 

> 
O 

(N 



Oh 



> 
OS 

oo 

00 



Elaine Goode 

University of Southampton, School of Physics and Astronomy, Highfield, Southampton, S017 

IB J, United Kingdom 

E-mail: ejg4g08@soton.ac.uk 

Matthew Lightman 

Department of Physics, Columbia University, New York, NY 10027 USA 
Department of Physics, Washington University, St. Louis, MO 63130 USA 
E-mail: mlightman@physics . wustl . edu 

The A/ = 3/2A'-> kk decay amplitude is calculated on RBC/UKQCD 32 3 x 64, L s = 32 dy- 
namical lattices with 2 + 1 flavors of domain wall fermions using the Dislocation Suppressing 
Determinant Ratio and Iwasaki gauge action. The calculation is performed close to the physi- 
cal pion mass (m n = 142.9(1.1) MeV) and with a single lattice spacing (cr 1 = 1.375(9) GeV.) 
We find Re(A 2 ) = (1.436 ± 0.063 stat ± 0.258 syst ) x l(T 8 GeV and Im(A 2 ) = (-6.29±0.46 stat ± 
1.20 syst ) x 10~ 13 GeV. These results are combined with the experimental result for e'/e to pre- 
dict Im(Ao) = -5.32(64) sta t(71) S y St x lCT n GeV within the Standard Model. We also perform 
a reweighting analysis to invesitgate the effects of partial quenching in the light-quark sector 
of our calculation. Following reweighting we find Re(A 2 ) = (1.52 zb 0. 14 s t a t) x l(T 8 GeV and 
Im(A 2 ) = (— 6.47±0.55 s tat) x 10~ 13 GeV, which are consistent with our main results. 
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1. Introduction 

The calculation of K — > %% decay amplitudes is motivated by a desire to understand the 
AI = 1/2 rule and CP violation in kaon decays. Such a calculation is a non-perturbative prob- 
lem requiring lattice techniques to make progress. Previous lattice calculations have relied on the 
quenched approximation and uncertain chiral extrapolations [1, 2, 3, 4]. In this talk we present the 
results from the first realistic lattice calculation of a K — > tctc decay amplitude, where we simulate 
the two-body decay directly on the lattice at nearly-physical kinematics. 

We proceed by evaluating matrix elements of the AS = 1 effective Hamiltonian 

#eff = Ow£Q& (1-1) 

i 

where Q are Wilson coefficients and <2; are four-quark operators. In this talk we consider only 
the AI = 3/2 transition, in which case only three operators contribute in equation (1.1). We find it 
convenient to evaluate unphysical K + — > 7i + 7i + matrix elements of the following operators: 

2(27,1) = is'd^d^L, {2( M ) = {jrf) L &d 3 ) R and 2 (8 , 8)mx = {fd^uU 1 )^ (1.2) 

where the operators are labelled according to their transformation under SU (3)l x SU (3)r, and the 
labels i and j on the quark fields label colour. These are related to the physical K + — > tc + tc° matrix 
elements via the Wigner-Eckart theorem. 

2. Details of the Simulation 

The analysis is performed on a single ensemble of 2 + 1 flavour domain wall fermions (DWF) 
with Dislocation Suppressing Determinant Ratio (DSDR)+Iwasaki gauge action at j8 = 1.75. The 
lattice size is 32 x 64 and the extent of the fifth dimension is L s = 32. The inverse lattice spacing 
is a -1 = 1.375(9) GeV. The ensemble is genereated with sea-quark masses am/ = 0.001 and am/, = 
0.045, corresponding to a unitary pion mass of approximately 170 MeV. We find the residual mass 
to be am ies = 0.00184(1). The correlation functions are calculated with valence quark masses of 
mi = 0.0001 and m s = 0.049, corresponding to a valence pion mass of m n = 142.9(1.1) MeV 
and kaon mass of yhk = 511.3(3.9) MeV. A total of 63 gauge configurations, each separated by 8 
molecular dynamics time units, are included in the analysis. 

Quark propagators with periodic and antiperiodic boundary conditions in the time direction 
were computed on each configuration with a source at t = 0. They were then combined so as to 
effectively double the time extent of the lattice. Meson correlation functions formed using the av- 
erage of the propagators with periodic and antiperiodic boundary conditions can be interpreted as 
containing forward propagating mesons originating at time t = 0, whereas those calculated with the 
antisymmetric combination can be interpreted as containing backward propagating mesons origi- 
nating from a source at t = 64. Strange-quark propagators, also with Periodic + Antiperiodic 
combinations, were generated with sources at tx = 20, 24, 28, 32, 36, 40 and 44 in order to cal- 
culate K — )• Tin correlation functions with kaon sources at these times, while the two-pion sources 
remained at either t = or t = 64. Thus we could achieve time separations between the kaon and 
two pions of 20, 24, 28 and 32 lattice time units in two different ways which increased the statistics. 
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These separations were chosen so that the signals from the kaon and two pions did not decay into 
noise before reaching the four-quark operator £?,. 

For physical K — > %% decays in the CM frame, the final state pions have equal and opposite 
non-zero momentum. We achieve this for the n + 7i + final state by giving momentum to the d-quark. 
The u- and s-quark propagators are generated with Coulomb gauge fixed wall sources and periodic 
spatial boundary conditions. Similarly, the d-quark propagators used in the zero-momentum pion 
and two-pion correlation functions are computed with Coulomb gauge fixed wall sources and pe- 
riodic spatial boundary conditions. However, the d-quark propagators used in the two-pion and 
correlation functions with non-zero momentum pions are generated with antiperiodic spa- 
tial boundary conditions and cosine sources. By imposing antiperiodic spatial boundary conditions 
the allowed quark momenta are p n = (71 + 2nn) /L, where L is the spatial extent of the lattice, cor- 
responding to a ground-state momentum of %jh in the direction in which antiperiodic boundary 
conditions have been used. In practice we impose antiperiodic boundary conditions in two spatial 
directions, which allows us to simulate pions with ground state momentum ±y/2n/L. This deci- 
sion is motivated by the expectation that pions with p = \/2n/L will correspond to a two-pion final 
state with energy close to mx ■ The use of cosine sources in K + — > 7i + n + decays is described in [5]. 

3. Analysis 

We extract the K — > %% matrix element j$ by fitting a constant to the left hand side of (3.1) 



C' Kjt7l is the K — > Tin correlator with a kaon source at t& , i labels the four-quark operator Q, which 
is inserted at time t, and Zk and Z nK are calculated from the kaon and two-pion correlators respec- 
tively, whose sources are at t = 0. The left hand side of equation (3.1) is plotted in Figure 1 for 
each operator. The figure demonstrates that sufficiently far from the kaon and two-pion sources we 
are justified in fitting to a constant. The fit results for (ZkZ kk ) are indicated on the plot. 

We also use a quotient method to extract the two-pion energy, as we find this improves the 
statistical precision of the fits. We fit the quotient of correlators C nn /(C n ) 2 ~ Ae _A£f to extract 
AE = E %% — 1E % . We then get the two-pion energy by calculating AE + 2E n , where in the case of 
p = 0, E n = m n while for p = y/ln/L, E n is found from a 2 parameter fit to the pion correlation 
function which also has p = \fj.%jh. The numerical results for all the meson masses and energies 
which we extract from the correlation functions are given in Table 1. From Table 1 we see that the 
kaon mass is not exactly equal to the two-pion energy, so our calculation is not quite on-shell. This 
is taken into account when estimating the systematic error on the final results. 

The finite volume matrix elements are related to the infinite volume amplitudes ^ (where i 
labels the 4-quark operator) using the Lellouch-Luscher factor [6, 7] . In particular we have 



where the quantity in square brackets contains the effects of the Lellouch-Luscher factor beyond 
the free field normalization, 8 is the s-wave phase shift, q n is a dimensionless quantity related to 



(3.1) 



C K (t K -t)C nn (t) Z K Z, 




(3.2) 
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units 


m% 


m K 


E„ t 2 




Exn,2 


hik — E n n,2 


lattice 
MeV 


0.10395(32) 
142.9(1.1) 


0.37193(91) 
511.3(3.9) 


0.1737(14) 
238.8(2.4) 


0.20948(63) 
288.0(2.2) 


0.3583(35) 
492.6(5.5) 


0.0136(35) 
18.7(4.8) 



Table 1: Results for meson masses and energies. The subscripts 0, 2 denote p = 0, p = \ 2n/L respectively. 




(a) (27, 1) operator (b) (8, 8) operator (c) (8, 8) mix operator 

Figure 1: K—> iziz quotient plots for p = \^2n/L. The two pion source is at t = while the kaon source is 
at t = 24. The dashed line shows the error on the fit. 

the individual pion momentum k n via q n = k n Lj2% and is a kinematic function defined in [6]. 
The pion momentum k n is calculated using the dispersion relation E nn = 2 ^/k^ + in^., and differs 
from y/2n/L due to interactions between the two pions. Once q n is known, 8 can be calculated 
using the Liischer quantisation condition [8], nn = 8(k K ) + <§>{qn)- As discussed in [5], the phase 
shift derivative is calculated using the phenomenological curve of [9]. This is necessary because 
we only have two values of the two-pion energy from which to extract the phase shift. 
The amplitudes ^ are related to the physical decay amplitude Ai via 

A 2 = a- 3 J~G F V ud V*Y l GW z ijW*h ( 3 - 3 ) 

V Z U 

where C; are the Wilson coefficients and Z, 7 are the renormalization constants. The Wilson co- 
efficients must be evaluated at the same scale and scheme as the renormalization constants. The 
renormalization constants are first evaluated in the RI-SMOM(^,^) scheme [10]. In order to min- 
imize discretization effects, this procedure takes place at a relatively low energy Hq = 1.145 GeV. 
A non-perturbative step-scaling function is then used to convert these results to a scale of 3 GeV, at 
which point a perturbative matching to the MS-NDR scheme is possible. The Wilson coefficients 
are known in the NDR scheme at the W-mass scale, and can be perturbatively run to the desired 
matching point of 3 GeV [11]. 

4. Results 

We calculate A 2 for the four different separations between the kaon source and two-pion 
source. Our final result, presented in equation (4.1), is an error weighted average over these four 
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results. The first error in equation (4. 1) is a statistical error, where the statistical uncertainties in the 
amplitude A, and lattice spacing (4% in total) are combined in quadrature with the statistical error 
on the renormalization constants (0.8% for ReA2 and 6% for Im(A2)). 

Re(A 2 ) = (1.436±0.063 st at±0.258 sys t) x 10~ 8 GeV, 

(4-1) 

Im(A 2 ) = (-6.29±0.46 stat ± 1.20 syst ) x 10~ 13 GeV. 

The second error in equation (4. 1) is systematic. The systematic errors in our calculation, for 
(Re(A2), Im(A2)) respectively, are from lattice artifacts (15%, 15%), uncertainty in the phase shift 
derivative (0.32% 0.32%), finite volume (6.2% and 6.8%), partial quenching (3.5% and 1.7%), 
uncertainties in the renormalization procedure (1.7% and 4.7%), unphysical kinematics (3.0% and 
0.22%), and perturbative truncation in the evaluation of the Wilson coefficients (7. 1 % 8. 1 %). Com- 
bining in quadrature, we find the systematic errors to be 18% for Re(A2) and 19% for Im(A2). 
Further details on how these errors are estimanted can be found in [12, 13]. 

5. Reweighting 

We use the technique of reweighting to test the consequences of the partial quenching in the 
light-quark sector of our calculation. The reweighting is performed in 30 increments from the 
simulated mass mf^ = 0.001 down to a value of mf a = 0.0001, corresponding to the valence light- 
quark mass. The results are shown in Figure 2. The rightmost point in Figure 2(a) shows the result 
for Re(A 2 ) before reweighting, while the remaining points show the results after reweighting to the 
mass indicated on the *-axis, ending with mf a = 0.0001 for the leftmost point. Similarly Figure 
2(b) shows the effects of reweighting on Im(A 2 ). Examining the figures, it can be seen that the 
errors on Re(A2) and Im(A2) grow, but the central values remain unchanged within the errors. Since 
reweighting effectively reduces the number of configurations contributing to the observables [14] it 
is natural that the statistical error should increase. However, the observation that the central values 
are unchanged confirms that partial quenching in the light quark does not introduce a significant 
source of systematic error. 

The final results after reweighting are shown in Table 2 where they are compared with the 
results before reweighting. 





mi = 0.001 


mi = 0.0001 (reweighted) 


Re(A 2 ) 
Im(A 2 ) 


1.436(63) x 10- 8 GeV 
-6.29(46) x 10~ 13 GeV 


1.52(14) x 10~ 8 GeV 
-6.47(55) x 10~ 13 GeV 



Table 2: A2 before and after reweighting. 



6. Prediction for Im(Ao) 

Assuming isospin symmetry, the CP- violating parameter e'/e can be expressed in terms of 
Re(Ao), Im(Ao), Re(A2) and Im(A2) according to equation (6.1): 



Re 



CO 



y/2\e\ 



Im(A 2 ) Im(A Q ) 
Re(A 2 ) Re(Ao) 



(6.1) 
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1.4 
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00001 (rw) 



0005 
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0001 (rv 



0005 
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0.001 (sin 



(a) Reweighting Re(A2) (b) Reweighting Im(A2) 

Figure 2: Reweighting A 2 from mf* = 0.001 to mf a = 0.0001. 



where ft) = Re(A2)/Re(Ao). Re yjj> l £ l> an d Re(Ao) are known experimentally and presented 
in Table 3. Combining these known factors with our lattice result for Im(A2)/Re(A2) we can 
determine the the unknown quantity Im(Ao) within the Standard Model, finding 



Im(A ) = -5.32(64) stat (71) syst x KT 11 GeV 



(6.2) 



The error on ImAo is obtained by combining the errors on the quantities in Table 3 in quadrature. In 
equation (6.3) below we compare the relative contribution to Im(Ao) /Re(Ao) from Im(A2) /Re(A2) 
and the term containing the experimentally known contributions: 



Im(A ) 
Re(A ) 

-1.60(19) stat (21) syst x lfT 4 



Im(A 2 ) 
Re(A 2 ) 

-4.38(34) stat (95) syst x lfT 5 



y/2\e\ e' 



CO 



(6.3) 



1.16(18) x 10 



-4 



Thus we see that while the error on the determination of Im(Ao) is dominated by the uncertainty 
in the experimental value of e'/e, the contribution of Im(A2)/Re(A2) to Im(Ao) is significant (about 
25% in the determination of Im(Ao) /Re(Ao)). 



Re(e'/e) 


(1.65 ±0.26) x 10~ 3 


ft) 


0.04454(12) 


|e 


(2.228 ±0.011) x 10" 3 


Re(A ) 


3.3201(18) x 10~ 7 GeV 


Im(A 2 )/Re(A 2 ) (lattice) 


-4.38(34) stat (95) syst x lO" 5 



Table 3: Experimental values of the components of equation (6.1) used in the determination of Im(Ao), 
together with the results for Im(A2)/Re(A2) from these proceedings. 
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7. Conclusions 

We have presented preliminary results for the Al = 3/2 K — > %% decay amplitude on 32 
lattices with 2+1 flavours of DWF and the Iwasaki-DSDR gauge action. We find 

Re(A 2 ) = (1.436±0.063 st at±0.258 syst ) x 1(T 8 GeV, 
Im(A 2 ) = (-6.29±0.46 stat ± 1.20 syst ) x 10~ 13 GeV. 

Our result for Re(A2) is in good agreement with the experimental result of 1.479(4) x 10~ 8 GeV 
obtained from K + decays. In the future we plan to undertake similar calculations of Re(Ao) and 
Im(Ao) [15], allowing e'/e to be calculated from first principles for the first time. In the mean 
time we make the prediction, based on the Standard Model, that Im(Ao) = — 5.32(64) sta t(71) S ys x 
10 _ 11 GeV. 
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